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Abstract 

The nonperturbative nature of inter-nucleon interactions is explored by varying 
the momentum cutoff of a two-nucleon potential. Conventional force models, which 
have large cutoffs, are nonperturbative because of strong short-range repulsion, the 
iterated tensor interaction, and the presence of bound or nearly-bound states. But 
for low-momentum interactions with cutoffs around 2 fin -1 , the softened potential 
combined with Pauli blocking leads to corrections in nuclear matter in the particle- 
particle channel that are well converged at second order in the potential, suggesting 
that perturbation theory can be used in place of Brueckner resummations. Calcu- 
lations of nuclear matter using the low- momentum two-nucleon force Vjowfc with a 
corresponding leading-order three-nucleon (3N) force from chiral effective field the- 
ory (EFT) exhibit nuclear binding in the Hartree-Fock approximation, and become 
less cutoff dependent with the inclusion of the dominant second-order contributions. 
The role of the 3N force is essential to obtain saturation, and the contribution to 
the total potential energy is compatible with EFT power-counting estimates. 



1 Introduction 



Conventional wisdom among nuclear physicists, as summarized by Bethe in 
his review of over 30 years ago [1] , holds that successful nuclear matter cal- 
culations must be nonperturbative in the inter-nucleon interactions. The pos- 
sibility of a soft potential providing a perturbative solution to the nuclear 
matter problem was discarded at that time, and saturation firmly identified 
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with the density dependence due to the tensor force [1]. Subsequent work on 
the nuclear matter problem [2] has not significantly altered the general per- 
spective or conclusions of Bethe's review (although the role of three-nucleon 
(3N) forces has been increasingly emphasized). However, recent results from 
renormalization-group-based low-momentum potentials, coupled with insight 
from effective field theory (EFT), mandate that these conclusions be revis- 
ited. These developments affect not only the technical conclusions but also 
the underlying philosophy of nuclear matter calculations. 

Nonperturbative behavior in the particle-particle channel for nuclear forces 
arises from several sources. First is a strongly repulsive short-range interaction, 
which requires at least a summation of particle-particle ladder diagrams [1]. 
Second is the tensor force, e.g., from pion exchange, which is highly singular 
at short distances, and requires iteration in the triplet channels [3,4]. Third 
is the presence of low-energy bound states or nearly-bound states, which are 
found in the nucleon-nucleon (NN) S-waves. These states imply poles in the 
scattering T matrix that render the perturbative Born series divergent. All 
of these nonperturbative features are present in conventional high-precision 
NN potentials, such as the Argonne t> 18 potential [5], which has been used 
in the most accurate ab initio calculations of nuclei and nuclear matter to 
date [6,7,8]. 

The philosophy behind the standard approach to nuclear matter is to attack 
these features head-on. This attitude was succinctly stated by Bethe [1]: 

"The theory must be such that it can deal with any NN force, including 
hard or 'soft' core, tensor forces, and other complications. It ought not to 
be necessary to tailor the NN force for the sake of making the computation 
of nuclear matter (or finite nuclei) easier, but the force should be chosen on 
the basis of NN experiments (and possibly subsidiary experimental evidence, 
like the binding energy of H 3 )." 

In contrast, the EFT perspective stresses that the potential is not an observ- 
able to be fixed from experiment (there is no "true potential"), but that an 
infinite number of potentials are capable of accurately describing low-energy 
physics [9] . In order to be predictive and systematic, an organization ( "power 
counting" ) must be present to permit a finite truncation of possible terms in 
the potential. If a complete inter-nucleon potential is used, including many- 
nucleon interactions, then all observable results should be equivalent up to 
truncation errors. The EFT philosophy applied to nuclear matter implies using 
this freedom to pick a convenient and efficient potential under the conditions 
of interest. 

A particularly useful class of energy- independent NN potentials is character- 
ized by a momentum cutoff A, which limits the resolution of details in the 
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potential. One starts with an accurate NN potential (such as Argonne v±s [5] 
or a chiral N 3 LO potential [10,11]) and a sufficiently high "bare" cutoff (all 
potentials are cut off at high momentum by form factors or other regulators) 
and evolves the cutoff to lower values, while requiring that observables (e.g., 
phase shifts) for external momenta up to the cutoff are unchanged [12,13]. 
This can be accomplished by renormalization-group (RG) equations that pre- 
serve the half-on-shell T matrix elements (and therefore the physical on-shell 
amplitude) or, equivalently, by Lee-Suzuki transformations [12,13,14] (effec- 
tive interactions in momentum space have also been derived in Ref. [15]). 
Varying the cutoff can be used as a powerful tool to study the underlying 
physics scales, to evaluate the completeness of approximate calculations, and 
to estimate truncation errors from omitted higher-order contributions. 

These variable-cutoff potentials reveal the resolution or scale dependence of 
the first two sources of nonperturbative behavior, which are dampened as 
high-momentum intermediate states are eliminated. In free space, the third 
source of nonperturbative behavior remains independent of the cutoff because 
the pole positions of shallow bound states that necessitate fine tuning are 
physical observables. However, this fine tuning is eliminated in the medium at 
sufficiently high density. Our results show that a low-momentum two-nucleon 
force ("Viowfc") eliminates all three sources of nonperturbative physics for bulk 
properties of nuclear matter. In short, a repulsive core is not constrained by 
phase shifts and is essentially removed by even a moderately low-momentum 
cutoff, the tensor force is tamed by a sufficiently low cutoff, and the shallow 
bound states become perturbative as a result of Pauli blocking. For cutoffs 
around 2fm _1 , which preserve phase shifts up to 330 MeV laboratory energy, 
the Born series in nuclear matter is well converged at second order in the po- 
tential, suggesting that perturbation theory can be used in place of Brueckner 
resummations. 

While evolving a soft potential from higher momentum is a new development in 
nuclear physics [12,13], attempts to use soft potentials for nuclear matter were 
common in the 1960's and early 1970's [16]. It had long been observed that 
a strongly repulsive core is not resolved until eight times nuclear saturation 
density [1]. Thus, saturation is not driven by a hard core (unlike liquid 3 He). 
However, these soft potentials were abandoned because they seemed incapable 
of quantitatively reproducing nuclear matter properties. Their requiem was 
given by Bethe [1]: 

"Very soft potentials must be excluded because they do not give saturation; 
they give too much binding and too high density. In particular, a substantial 
tensor force is required." 

From the EFT perspective, a failure to reproduce nuclear matter observables 
should not be interpreted as showing that the low-energy potential is wrong, 
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but that it is incomplete (e.g., many-body forces). This misconception still per- 
sists and has recently led to the conclusion that low-momentum two-nucleon 
interactions are "wrong" since they do not give saturation in nuclear matter 
and finite nuclei are overbound for lower cutoffs [17,18]. The missing physics 
that invalidates this conclusion is many-body forces, which were completely 
neglected in those studies. 

In a low-energy effective theory, many-body forces are inevitable; the relevant 
question is how large they are. It is established beyond doubt that for all 
realistic potentials, a significant three-body force is required to describe light 
nuclei [6,7,19,20]. For variable-cutoff potentials, the three-body (and higher 
many-body) components of the potential also evolve with the resolution scale. 
A full RG evolution of the combined two- and three-body potential is not 
yet available, but can be approximated by fitting the three-body potential at 
each cutoff to the form of the leading-order 3N force from chiral EFT. The 
first study of this sort revealed that the leading three-body components of 
the 3N force become perturbative at lower cutoffs [21], which implies they are 
tractable in many-body calculations. In this paper, we apply these potentials 
to nuclear matter and study the convergence properties of low-momentum 
interactions. 

The plan of this paper is as follows. In Section 2, a variable-cutoff potential is 
applied to the Born series for the two-nucleon amplitude in free space and as a 
function of density. The convergence properties as a function of cutoff are an- 
alyzed quantitatively using an approach introduced long ago by Weinberg. In 
Section 3, calculations of nuclear matter using the low-momentum two-nucleon 
interaction V\ ow k with a corresponding leading-order 3N force from chiral EFT 
are presented in the Hartree-Fock approximation and an approximation to sec- 
ond order. Nuclear saturation naturally arises at the Hartree-Fock level and 
becomes less cutoff dependent with the inclusion of the dominant second-order 
contributions. Our conclusions are summarized in Section 4 along with a plan 
for further calculations of nuclear matter. 



2 Perturbative ladders with low-momentum interactions 

In this section, we demonstrate that two sources of nonperturbative behavior 
in the particle-particle channel, namely "hard core" 1 scattering and iterated 
tensor force contributions, can be rendered perturbative by using the RG to 
lower the momentum cutoff. We also show that, at sufficiently high density, 
Pauli blocking suppresses the nonperturbative physics due to bound or nearly- 

1 We use "hard core" as a shorthand for the strong (but not infinite) short-range 
repulsion in conventional nuclear forces. 
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bound states [22,23]. In the language of EFT, this means that the infrared- 
enhancements that occur in the presence of such states in free space [24] are 
absent in nuclear matter. 

In order to keep the presentation self-contained, we briefly review the RG 
method used to construct low-momentum NN interactions ("Vi ow fc") starting 
from an arbitrary potential model Vnn [12,13,14]. In a given partial wave, 
two-nucleon scattering is described by the T matrix, which with standing- 
wave boundary conditions and in our conventions (h = ton = 1) is given by 2 

r(V u u*\ t/ (V »x 2 ^r 2j VNN(k',q)T(q,k;k 2 ) 

T(k ,k;k ) = Vn N (« ,k) + -V q dq — . (1) 

n Jo fc — q 

Here V denotes a principal value integration and the scattering phase shifts are 
given in terms of the diagonal T matrix elements, tan8(k) = —kT(k,k;k 2 ). 
Next, we define a low-momentum version of Eq. (1) in which a cutoff A is 
imposed on loop integrals and the external legs of the interaction, 

rp ci i i2\ t t ci i\ . ^ /" A 2 i Viowk(k', q) T\ ow k(q, k; k 2 ) 

T lowk (k ,k;k ) = Vi owk {k ,k) + -V q dq — . (2) 

n Jo k — q l 

As the matching condition, we demand that the half-on-shell T matrices are 
equivalent T\ ow j,(k' , k; k 2 ) = T(k' , k; k 2 ) for k',k < A and all cutoffs. This 
ensures that V\ ow k gives the same low- momentum two-nucleon observables as 
Vnn an d has the advantage that it avoids energy- dependent interactions, which 
are inconvenient in many-body applications. The cutoff independence of the 
low-momentum T matrix leads to an RG equation (for details see [14]) 

d T/ M 2 \/ lowfe (A; , ,A)T lowfe (A,A;;A 2 ) 
IK V ^ k{k > k) = n 1 - (k/A) 2 • (3) 

Whereas the T matrices in Eqs. (1) and (2) are right-side half-on-shell, the T 
matrix in the RG equation is left-side half-on-shell and is not RG invariant. 
This RG equation is asymmetric in k and k', so it generates a non-hermitian 
Viowfc- We avoid this by working with a symmetrized version of Eq. (3), 

d_, r ,i, i, _ 1 ( Vi owk (k',A)T lowk (A,k;A 2 ) T lowk {k', A; A 2 ) V lowk {A, k) \ 
dA V, owk {K,K)- ^ l_(jfc/ A )2 + l-(jfe'/A)2 )' 

(4) 

which preserves the on-shell T matrix (i.e., observables) and is equivalent to 
the Okubo hermitization of effective interaction theory. We note that Vi ow fc 
can be obtained by numerically integrating the RG equation with Vnn as the 
large-cutoff initial condition, or equivalently by Feshbach projection-operator 
techniques such as the Lee-Suzuki method (for details see [13,14]). 



2 This quantity, which is real, is sometimes referred to as the K matrix. 
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2.1 Born series in free space vs. in medium 



We can now use the RG to study the convergence properties of free-space and 
in-medium T matrices as we change the resolution scale A. In perturbation 
theory, the T matrices are given by the Born series, 

T(E) = V + VG (E)V + VG (E)VG (E)V + ..., (5) 

where the unperturbed propagator G (E) is 

G (E) = 1 (free-space) , G (E) = ^ F (in-medium) , (6) 
hi — sIq hi — Hq 

and Qk F is the Pauli-blo eking operator that allows scattering only to unoccu- 
pied states above the Fermi momentum hp. In momentum space, the first and 
second Born terms are given by 



T {1 \k',k;E) = V(k', k) , 

T< 2 > (k', k;E) = -V I q 2 dq V ^ V{ ?> k) (free-space) , (7) 
Ti J E — q l 

T^,k;E)^-V /rt ^f^jy* 1 (in-medto), 

where P is the total pair momentum, and for S-waves the Pauli-blocking 
operator is given after angular integration by 



Q(q,P;k F ) = (q,P\Q kF \q,P) 



for q < yfj% - P 2 /4 

1 for q > k F + P/2 (8) 
q + P /4 k_p ^j ierw j se _ 



qP 



In the above equations, V(k', k) stands for either Vi ow fc or the input potential 
Vnnj an d the integration is from to A or to oo, respectively. For simplicity, 
we neglect self-energy effects and restrict our attention to total pair momentum 
P = in this section. The dependence on P is generally very weak, and so 
our conclusions hold for non-zero P as well. 

The left panel of Fig. 1 shows the first- and second-order contributions to the 
diagonal T matrix elements calculated from the Argonne i>i 8 potential in the 
1 S channel. It is evident that perturbative expansions of the free-space and 
in-medium T matrices are not possible, as the second-order contribution is sev- 
eral times larger than the first-order contribution over all momenta considered. 
This divergence continues in higher orders. The strong repulsive core in the 
Argonne v is interaction, or any other conventional NN potential model, leads 
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Fig. 1. Perturbative contributions to the diagonal So T matrices calculated with 
the Argonne v\% potential (left panel) and with Vj Q wfc f° r A = 2.1fm _1 (right 
panel). The first-order interactions are given by the solid lines. The second-order 
free-space (in- medium) contributions are denoted by the dashed (dot-dashed) lines. 
The in-medium contributions are evaluated at saturation density and for total pair 
momentum P = 0. Here and in the following Vjowfc is calculated from the Argonne 
vis potential. 

to strong high-momentum matrix elements, and the momentum-independent 
nonperturbative behavior observed in Fig. 1 is a signature of "hard core" 
scattering. We emphasize that cores are not constrained by fits to low-energy 
two-nucleon data, and therefore "hard core" divergences of perturbation the- 
ory are not dictated by physics. 

Pauli blocking does not change this picture, and in nuclear matter the second- 
order contribution is comparable to the free-space one, as can be seen from 
Fig. 1. Therefore, "hard core" scattering always dominates. One might have 
expected that the second-order contribution is small at threshold in nuclear 
matter, since Pauli blocking suppresses the infrared enhancements that occur 
in free space at low energies. Clearly, this is not the case with conventional 
potential models due to the cores. 

Referring to the right panel of Fig. 1, we find that the results obtained from 
Mowfc with A = 2.1 fm -1 exhibit a very different pattern. We first observe 
that the second-order contributions to the amplitudes are smaller than the 
first-order Born term over the range of momenta considered. The second- 
order free-space contribution is largest at zero energy due to the large 1 So 
scattering length, 3 while it decreases rapidly away from threshold, indicating 
that "hard core" scattering is absent. We emphasize that the diagonal free- 

3 With the normalization used here, the contributions at k = will sum to the 1 So 
neutron-proton scattering length of ai So = — 23.7fm. 
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A: \P/2±k\>k l7 and Ikl < A 

F 




Fig. 2. Pauli-allowed intermediate-state phase space (hatched region "A") for scat- 
tering in the particle-particle channel. Nucleons are initially inside the Fermi sea 
(shaded region "F" ) and interactions are restricted to momenta below the cutoff A. 

space T matrix elements obtained from an exact solution using either the 
Argonne v\% potential or V\ ow k are identical up to 330 MeV lab energy. This 
demonstrates that the convergence behavior of the Born series in free space is 
strongly scale-dependent, with lower cutoffs being advantageous. 

For the in-medium calculation, we find the striking result that the second-order 
contribution is very small over all momenta. This dramatic improvement of 
in-medium convergence can be easily understood from a simple phase-space 
argument. In the present calculation for total pair momentum P = 0, the 
limits of the intermediate-state integrals in Eq. (7) are from k F to A for V\ ow k- 
In contrast, the loop integrals for conventional potentials require large upper 
limits of 15-20 fm _1 . Therefore, the restriction of phase-space substantially 
reduces the size of higher-order corrections. This reduction is also aided by 
weaker interactions at higher momenta in the dominant S-waves. 

This phase-space argument also holds for general P ^ 0. In Fig. 2, we show 
the allowed phase space for the general case, where a pair of particles with 
total momentum P initially in the Fermi sea (the shaded region "F") scatters 
to Pauli-allowed intermediate state (the hatched region "A"). It is clear from 
Fig. 2 that lowering the cutoff suppresses the intermediate-state phase space 
for P ^ case as well. When P is very close to 2k F the suppression is less 
extreme, but a corresponding Vy ow kGViowk curve as in Fig. 1 is still less than 
0.3 fm in magnitude throughout the range of k. 

However, the cutoff also has to be greater than the largest relative momentum 
with non-negligible contributions to the observable of interest. It is evident 
from Fig. 2 that for large total pair momentum (e.g., P = 2kp), cutoffs A < 2/cp 
cut into the Fermi spheres. For general interactions, taking the cutoff too small 
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could complicate the task of maintaining cutoff-independent observables. This 
is because the reshuffling of contributions from loop integrals to the effective 
interaction in the RG equation is defined in free space. In this paper we take a 
pragmatic approach to this issue. Energy-per-particle results for A < 2k F are 
found to be very weakly dependent on the cutoff, which indicates that such 
configurations contribute negligibly to bulk properties. 4 Moreover, region "F" 
vanishes as P — > 2/cp and the average pair momentum in the Fermi sea is 
-Pav — \j6/5kp. Therefore, it is likely that the limit A ^ 2k-p set by the 
P = 2kp configurations is too conservative. 

In summary, we find that the perturbative convergence of two-nucleon scat- 
tering is improved dramatically when we use the RG to evolve the interaction 
down to lower cutoffs. The resulting Vi ow k is by construction phase-shift equiv- 
alent to the input potential, but Vj ow fc removes the difficulties due to repulsive 
cores. In-medium ladders become perturbative because lower cutoffs restrict 
the intermediate-state phase-space integrations, and because NN scattering in 
the dominant S-waves is weaker at higher momenta. 

2.2 Weinberg eigenvalue analysis 

The preceding conclusions can be made rigorous by adopting a method intro- 
duced by Weinberg in Ref. [26] that provides quantitative conditions for the 
perturbative convergence of the T matrix. This allows us to identify the cut- 
offs at which "hard core" scattering and iterated tensor contributions become 
perturbative. Moreover, the analysis shows that Pauli-blocking removes non- 
perturbative behavior due to low-energy bound or nearly-bound states. Thus, 
we will establish as a main result of this paper that the in-medium T matrix 
is perturbative over a wide range of cutoffs and densities. 

Following Weinberg [26], we study the spectrum of the operator G (z)V for 
complex energies z 

G (z)V\* v (z)) =ri v (z)\* v (z)), (9) 

where the index v labels the discrete eigenvalues and eigenvectors. As demon- 
strated in [26], the perturbative Born series for T(z) diverges if and only if 
there is an eigenvalue with ^(z)! ^ 1. The necessity of this condition is clear, 
since 

oo oo 

E v(G (z)v) n \Hf v (z)) = £ (ifc(s)) V \Hf v (z)) , (10) 

n=0 n=0 

4 We also note that observables for which relative momenta k > 2fm~ 1 contribute 
significantly will be model dependent. This is the case, for example, for P-wave 
pairing gaps for neutron Fermi momenta kp > 2fm _1 [25]. 
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and thus 



T(z) \* v (z)) = (l + Vu(z) + r, v {zf + ...)V \*„(z)) 



(11) 



diverges for ^ 1. Furthermore, Weinberg also demonstrates that the 

rate of perturbative convergence is controlled by the largest ^(z)!, with 
smaller values implying faster convergence. 

A rearrangement of Eq. (9) gives a simple interpretation of the eigenvalues 
r\ v {z) in terms of the Schrodinger equation, 



The eigenvalue r\ v (z) can thus be viewed as an energy-dependent coupling 
constant that must divide V to produce a solution to the Schrodinger equation 
at energy z. If V supports a bound state at z = B < 0, then there is some 
v with rj u (B) = 1, which implies a divergence of the Born series for nearby 
energies. However, what matters for convergence at a given energy z is not 
simply the presence of nearby physical bound states, but rather the complete 
set of eigenstates that can be shifted to z when the interaction is divided by 
7) v (z). At negative energies, a purely attractive V gives positive r\ v {z) values, 
while a purely repulsive V gives negative values, as the sign of the interaction 
must be flipped to support a bound state. For this reason, we follow convention 
and refer to negative eigenvalues as repulsive and positive ones as attractive. 
In the case of conventional NN interactions, the repulsive core generates at 
least one large and negative eigenvalue that causes the Born series to diverge. 

We apply the spectral analysis to the coupled 3 Si- 3 D! channel and show in 
Fig. 3 the evolution of the two largest Weinberg eigenvalues evaluated at the 
deuteron pole z = Bd, as we lower the cutoff. From the free-space curves (solid 
lines), we observe that the repulsive eigenvalue associated with the core and 
the tensor force is initially large but decreases with cutoff, with \r] u (B,i)\ < 1 
for A < 3.5 fm -1 . In contrast, the physical divergence corresponding to the 
deuteron pole remains invariant, r] u (B ( i) = 1, as it should. This suggests that 
the nonperturbative behavior due to the strong core and strong tensor force 
is only present for large-cutoff interactions. 

The finite-density curves (dashed lines) show a similar dependence on cut- 
off, but now the nonperturbative attractive eigenvalue associated with the 
deuteron in free space is tamed by Pauli-blocking effects. This result is general, 
as illustrated by Fig. 4, where we show the density dependence of the Wein- 
berg eigenvalues for various cutoffs. As the density increases, Pauli-blocking 
drives the largest attractive Weinberg eigenvalue to a perturbative regime. 

It is important to note that the repulsive eigenvalue is significantly reduced 
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Fig. 3. Evolution of the two largest Weinberg eigenvalues with A in the coupled 
3 Si- 3 Di channel in free-space (solid) and at saturation density (dashed). The lighter 
curves correspond to attractive eigenvalues and the darker curves to repulsive eigen- 
values. 
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Fig. 4. Dependence on density of the two largest Weinberg eigenvalues in the 
coupled 3 Si~ 3 Di channel. 

when the cutoff is lowered from A = 3fm _1 to A = 2fm _1 . Second-order ten- 
sor contributions strongly excite intermediate-state momenta peaked about 
k ~ 2.5-3.0 fm -1 in nuclear matter [27]. Consequently, when the cutoff is low- 
ered below this range, the strength of the tensor force decreases, leading to a 
smaller Weinberg eigenvalue. Chiral EFT interactions have Weinberg eigen- 
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values similar to the A = 3 fm curves, and thus their in-medium convergence 
properties can be improved by using the RG to run their cutoff to lower values. 

We therefore reach a promising conclusion: Large-cutoff sources of nonpertur- 
bative behavior can be eliminated using the RG to lower the momentum cut- 
off, while the physical source of nonperturbative behavior due to bound and 
nearly-bound states is suppressed in the medium by the Pauli principle. For 
simplicity, the above analysis was for NN interactions only, but as we show 
in the next section, a consistent treatment of many-nucleon forces does not 
change the general conclusion. 



3 Nuclear matter results 

It is evident from the preceding section that for nuclear interactions with 
large cutoffs, the strong cores and tensor force necessitate a resummation of 
particle-particle scattering in nuclear matter. This resummation is typically 
carried out in the Brueckner approach, which sums iterated particle-particle 
scattering within a self-consistent mean- field (for reviews see [28,29]). 

In Fig. 5, we show the particle-particle contributions to the energy of sym- 
metric matter in perturbation theory. The calculations are carried out using 
a continuous single-particle spectrum with the effective mass approximation, 
and all intermediate-state phase-space integrations are angle-averaged. Fig. 5 
demonstrates that for large-cutoff interactions, such as the Argonne t>i 8 po- 
tential, iterated particle-particle scattering is nonperturbative. We repeat that 
this behavior is not required by nucleon-nucleon scattering data, and may be 
viewed as an artifact of the large cutoff in the Argonne v±s potential. In stark 
contrast, for the low-momentum interaction Vi ow fc with A = 2.1 far 1 we ob- 
serve that particle-particle ladders are perturbative in nuclear matter. These 
results are due to restricted phase space and follow from the behavior of the 
Weinberg eigenvalues, as discussed in the previous section. 

With two-nucleon low-momentum interactions only, particle-particle correla- 
tions do not lead to saturation within the density range where nuclear forces 
are well-constrained experimentally. We also note that for two-nucleon forces 
with cores, saturation occurs at Fermi momenta k-p > 1.5 fm" 1 in the Brueck- 
ner approach [28,29,30], in contrast to the empirical value of k F = 1.35 fm" 1 . 
Such calculations overbind nuclear matter by several MeV relative to the em- 
pirical value of E/A = —16 MeV. Moreover, it is well-established that 3N in- 
teractions present important contributions in light nuclei [7,19,20]. In Ref. [21], 
we showed that 3N contributions are smaller for low-momentum interactions, 
and that 3N forces are required to guarantee regulator-independent results. 
Therefore, it is imperative that 3N forces are included in calculations of nuclear 
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Fig. 5. Particle-particle contributions to the energy per particle in symmetric nu- 
clear matter for the Argonne v±s potential (thick lines) and the low-momentum 
interaction V\ ow k with A = 2.1 fin -1 (thin lines). 

matter. 

As discussed in the introduction, an RG evolution of the combined two- and 
three-nucleon interactions is not yet available. However, we have argued that 
this evolution can be approximated by fitting the leading-order 3N force from 
chiral EFT at each cutoff to be used with V\ ow k [21]. This is motivated by 
the fact that Viowfe becomes independent of the starting two-nucleon potential 
as the cutoff is lowered to A < 2fm _1 . The same mo del- independent V\ ow k is 
obtained by evolving a chiral N 3 LO two-nucleon interaction to smaller cut- 
offs. The RG evolution induces contact interactions which are of higher or- 
der in a chiral expansion and are necessary to maintain cutoff-independent 
two-nucleon observables. This indicates that Vi ow fc effectively parameterizes 
a chiral EFT interaction with higher-order contact terms, and supports our 
assumption that the low-momentum 3N interaction corresponding to VJ ow fc is 
well-approximated by the leading-order chiral 3N force. 

Therefore, we have augmented Vj ow fc with the leading-order chiral 3N force [21]. 
The latter contains a long-range 2-7r-exchange part V c , an intermediate-range 
l7r-exchange part Vd and a short-range contact interaction Ve [31,32]. The 
27r-exchange interaction is given by 

c 2\2fJ ^ k ( q ? + ml)( q ] + ml) ^ 




**.. AV18 - 
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where q« = k • — kj denotes the difference of initial and final nucleon momenta 
and 



7 a/3 
ijk 



4cim 2 2c_3 

f2 ' f2 

J TT J TT 



IT * • 



+ E T2 ^ T Z ■ x ' ( 14 ) 



f 2 



while the l7r-exchange and contact interactions are given respectively by 



V D = 

v E = 



9a c d 



E 



8 J 2 i^vfe Qj 

E ( T i ' T fe) 



{n ■ Tj) (<Ti ■ qj) , 



(15) 
(16) 



In applying Eqs. (13)— (16), we use qa = 1.29, f n = 92.4 MeV and m n = 
138.04 MeV and the q constants extracted by the Nijmegen group in a partial 
wave analysis with chiral 27r-exchange [33]. These are c\ = — 0.76 GeV -1 , 
c 3 = —4.78 GeV -1 and c 4 = 3.96 GeV -1 . The fit values for the c.d and c E 
low-energy constants are tabulated in Ref. [21] for A x = 700 MeV. 

Motivated by the above discussion, we proceed to calculate the energy per 
particle in symmetric nuclear matter from Viowfe and V3N in many-body per- 
turbation theory. The V\ ow k and V3N Hartree-Fock contributions are given by 



E 



(i) 

Mow; 



V 



Tr Tr 

J - L <T1,T1 - L1 (T2,T2 



J (2ti) 3 J 



dko 



n kl n k2 (12|y lowfe (l-P 12 )|12), 



(17) 



E 



(i) 
v 3N 



V 



Tr (Tl >Tl Tr CT2 )T2 Tr^ i 



' ra J (27p y (2^)3 y 



dk.2 f dk 3 



(2tt) 3 J (2tt) 



3 n fc2 n fc3 



x/ 2 (p,g)(123|\/ 3 N^i23|123), 



where the denote zero temperature Fermi-Dirac distributions, P\2 is the 
exchange operator for spin, isospin and momenta of nucleons 1 and 2, and V 
is the volume. Moreover, the momentum-conserving delta functions are not 
included in the NN and 3N matrix elements. For the 3N contribution, 



hip, <i) = ex P 



V + 3g 2 /4 N 
. A 2 . 



(19) 



is the regulator used in the 3N force fits in Ref. [21] and p and q are Jacobi 
momenta. Exchange terms are included by means of the antisymmetrizer 



-4-123 — (1 + P12P23 + ^13-^23) (1 — -P23) 

— 1 ~~ P12 ~ P13 ~ P23 + P12P23 + -P13-P23 



(20) 
(21) 



where the direct, single-exchange and double-exchange contributions are ap- 
parent in Eq. (21). The regulator h(p, q) is totally symmetric when expressed 
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in the nucleoli momenta kj, and thus all exchange terms contain the same 
regulator as the direct contribution. 



For symmetric nuclear matter, we obtain the following Hartree-Fock contri- 
butions from the leading-order chiral 3N force. For the contact part (Ve) and 
the l7r-exchange 3N interaction (V D ) we have 



E 



(i) 

E 



CE 



(2tt) 4 tfA x k F 



2k F y^'^T fc F + | 

J P 2 dP J p 2 dp J q' 2 dq' 



xF(p,P;k F )G(q',P; k F ) %(p,2q' /3) , 



(22) 



E 



(i) 

D 



9 g A c D 



2k f 



,2_P 2 



A (2tt)M/;* A X A$ 



J P 2 dP J p 2 dp J q' 2 dq' 



4p 2 



4p 2 + m 2 



xF(p,P;k F )G(q',P; k F ) flip, 2q'/3) , 



(23) 



where F and G are phase-space functions 



F(p,P;k F ) 



G(q',P; k F ) 



2 ^E- 



9 2 -P 2 /4 
~pP 



for p <k F - P/2 



otherwise 



k 2 -(q'-P/2) 2 

7p 



for q' < k F - P/2 
otherwise. 



(24) 



(25) 



It is useful to have simple estimates for these contributions. If the regula- 
tor is approximated by unity over the integration region, one has E^/A = 
— ce k F /(127r A A x ). This is reasonable because of the high power in the ex- 
ponential regulator /r. We further note that in the chiral limit, m n — > 0, the 
ratio of D- to _E-term contributions is E$ / 'E$ = —QaCd/^ce- This is an 
overestimate for finite pion mass by ~ m^jk\ « 25%. 
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1.6 
1.9 
2.1 
2.3 



CD 

2.080 
-1.225 
-2.062 



ce 
0.230 
-0.405 
-0.625 



-2.785 -0.822 

Table 1 ' 
Values for the low-energy couplings cd and eg adjusted to the 3 H and 4 He binding 
energies in [21]. As discussed in the text, the couplings for A = 2.1 fm -1 and A = 
2.3 fm -1 have been interpolated. 

The remaining 27r-exchange 3N interaction leads to 



E }) 



9 



2 2fe P 

U I P 2 dP 



pi 



k 2 - 

K F 4 



fc F + £ 



(2tt) 5 H 

r 

xT(p,P;k F )g(q',P; k F ) /|(p, 2q>/3) 
ki2 • k 2 3 



p 2 dp / q' 2 dq 



i 

d cos 9 P J d cos q / 

-l 



2tt 



X 



+ 



Cim 2 



+ 2. 



k 2 



ft 1(^2 + 0(^3 + (kh + ml)\ 

cs ( (k 12 ■ k 23 ) 2 fcj, \ 

2/4 ^ (A; 2 2 + ml)(P 23 + ml) (P 12 + ml) 2 J 



c 4 



k 12 x k 



23 ; 



2/ 4 V(^2 + 0(^3 + 



(26) 



where kjj = kj — k, and the terms with k 12 only are single-exchanges. Due 
to the angular dependence of the 27r-exchange 3N interaction, the angular 
integrations do not simplify and the phase-space restrictions are 



F(p,P;k F ) 



GW,P; fcp) 



1 

k 2 

| cos 9 P | < -E- 



1 

k 2 

COS^q/ < -E- 



p 2 - P 2 /4 
pP 



,'2 



q'P 



P 2 /4 



for p < k F — P/2 
otherwise 

for q' < k F - P/2 
otherwise. 



(27) 



f28) 



When the regulator is neglected, the In- and 27r-exchange 3N force contri- 
butions can be evaluated analytically or can be reduced to one-dimensional 
integrals. We give the corresponding expressions in the Appendix. 

We emphasize that the parameters for the 3N force were adjusted to the 3 H 
and 4 He binding energies in [21], and therefore the results presented in this 
section are predictions. To obtain results with A = 2.1 fm -1 and A = 2.3 fm -1 , 
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L 3N J 
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// 
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Fig. 6. Hartree-Fock (left figure) and Hartree-Fock plus dominant second-order con- 
tributions (right figure) calculated from V\ ow k and V3N for various cutoffs. Details 
of the approximate second-order calculations are given in the text. 

we have interpolated between the smoothly- varying c.d and ce fit values. 5 The 
values for the low-energy couplings cd and ce are listed in Table 1. Finally, 
we note that all Hartree contributions vanish for the leading-order chiral 3N 
interaction due to their spin-isospin structure. 

In Fig. 6, we present Hartree-Fock results for the energy per particle in sym- 
metric nuclear matter calculated from VJ ow fc and V3N for various cutoffs. The 
remaining integrals in Eqs. (22), (23), and (26) were calculated numerically. 
The Hartree-Fock results show that nuclear saturation is due to 3N forces 
when mo del- independent, low-momentum NN interactions are used. As can 
be seen from the Vj ow fc results in Fig. 6, the effects of 3N forces is very small 
at low densities k-p < 0.8 fm -1 . At the Hartree-Fock level, we find a minimum 
of E/A « —(2.5 — 8) MeV for Fermi momenta kp ~ (1.1 — 1.3) fm -1 over the 
cutoff range considered. It is quite promising that nuclear matter saturates in 
the Hartree-Fock approximation for low-momentum NN and 3N interactions, 
as conventional NN and 3N interactions require complicated nonperturbative 
treatments to achieve binding. However, we have not yet formulated a power 
counting appropriate for finite density that prescribes these ingredients at 
leading order. 

As can be seen from Table 2, the dominant 3N contributions are due to the 
27r-exchange interaction. For kp = 1.2 fm~ and A = 1.9 fm~ , the 3N Hartree- 



Fock expectation values are E$/A = 1.2 MeV, E^'/A = 
EW/A = 5.6 MeV. The repulsive single-exchange c-terms are more than a 
factor two larger than the attractive double-exchanges. The dominance of the 

5 For A = 2.5 fm -1 , we use the (b) fit of Ref. [21], since this lies smoothly between 



(i) 



-0.8 MeV and 



the A = 1.9 fm" 1 and A = 3.0 fm" 1 fits. 
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Hartree-Fock 



V e 



Hartree-Fock + dominant second order 
T V lowk V c V D V E 



1.0 



1.6 
1.9 
2.1 
2.3 



12.44 
12.44 
12.44 
12.44 



-19.62 
-18.18 
-17.35 
-16.56 



1.65 
1.67 
1.67 
1.67 



0.42 
-0.25 
-0.42 
-0.56 



-0.22 
0.40 
0.62 
0.81 



15.50 
16.29 
16.92 
17.60 



-26.58 
-26.81 
-27.04 
-27.27 



1.49 
0.85 
0.11 
-0.89 



0.34 
-0.09 
0.05 
0.43 



-0.29 
0.55 
0.79 
0.85 



1.2 



1.6 
1.9 
2.1 
2.3 




20.86 
21.80 
22.87 
24.32 



-37.66 
-37.38 
-37.53 
-37.95 



4.59 
3.99 
2.27 
-0.38 



1.03 
-0.50 
-0.37 

0.51 



-0.65 
1.28 
1.82 
1.78 



1.35 1.6 
1.9 
2.1 
2.3 



22.67 
22.67 
22.67 
22.67 



-42.47 
-38.82 
-36.74 
-34.77 



10.75 
11.95 
12.19 
12.30 



2.59 
-1.69 
-2.91 
-3.97 



-1.21 
2.34 
3.68 
4.89 



26.09 
26.75 
28.05 
30.06 



-47.85 
-46.72 
-46.47 
-46.45 



8.73 
9.14 
6.99 
3.10 



1.96 
-1.16 
-1.33 
-0.35 



-1.12 
2.24 
3.22 
3.26 



Table 2 

Expectation values of the kinetic energy (T), V\ owk and the different V3N contri- 
butions in MeV. The expectation values are obtained with the Feynman-Hellman 
method, and A and kp are given in fm _1 . 

single-exchange terms holds for other densities and cutoffs as well. 

Next, we compute the dominant second-order contributions to the energy per 
particle. The approximate second-order calculation is carried out in two steps. 
First, we convert the 3N force into a density-dependent NN interaction V^ 3 n 
by summing the third particle over occupied states in the Fermi sea, 

(12 I F 3N 1 1'2') = T W 3 J n fe3 / R (123) / R (l'2'3) (123 | ^ 3N I 1'2'3). (29) 

A first-order calculation of using V^n includes the single-exchange con- 
tributions of the Hartree-Fock calculation with Vw We then calculate the 
second-order contributions E^ from V\ owk P ms density-dependent K 3 n, 



E& 1 Jt / 

v - 4 11 I T W, 

i=i \ 



x 



(12 I (Vi 



low k 



(2tt) 3 / 
+ F 3 n)(1-Pi2)|34)| 2 



n kl n k2 (1 - n ks ) (1 - n k4 ) 



£ k3 + £fc 4 — e fci — Cfc 2 



(27r) 3 ^(k 1 + k 2 -k 3 -k 4 ). 

(30) 



For the intermediate-state integrations, the phase-space is angle-averaged and 
we use a continuous spectrum for e k = k 2 / (2m*). The angle-averaging approx- 
imation is expected to be reliable [34]. Here the effective mass is determined 
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kpffm ] 



Fig. 7. Ratio of V3N to Vi ow fc expectation values in the approximate second-order 
calculation. The expectation values are obtained by the Feynman-Hellman method. 
For nuclear matter we take for (V3N) the total 3N contribution, whereas the values 
shown for 3 H and 4 He are the largest 3N contributions (because of cancellations in 
the light nuclei). 



from the first-order self-energy correction from V\ ow k + ^3N at the Fermi sur- 
face, 

m* _ ( k dE{k) 
m \ 



m dk 



(31) 



k=kp 



where the spin-isospin independent part of the self-energy S(/c) is given by 



E(*) = \ Tr, 



Tr 

0"1,T1 J - J -CT2,T2 



dk 2 

W) 



n k2 (12 I (V lowk + V m ) (1 - P 12 ) I 12) . (32) 



In this approximation, we find for the effective mass m*/m = 0.72, 0.67 and 
0.65 for A = 1.9 fm -1 at kp = 1.0 fm~ , 1.2fm _1 and 1.35 fm _1 respectively. 



Our nuclear matter results including these dominant second-order contribu- 
tions are shown in the right panel of Fig. 6 and are promising in a number of as- 
pects. First, we observe that the second-order corrections move the saturation 
curves towards the empirical value. With the additional attraction, we find a 
minimum of E/A ~ —(11 — 12) MeV for Fermi momenta k-p w (1.2 — 1.35) fm _1 
over the cutoff range considered. Third-order particle-particle and hole-hole 
contributions are found to be small (< lMeV). Second, the cutoff depen- 
dence is dramatically reduced when second-order contributions are included. 
We note that the determination of the q are not unique [10,33,35,36]. For ex- 
ample, the values C3 = —3.2 GeV -1 and C4 = 5.4 GeV -1 determined by Entem 
and Machleidt [10] would yield an additional ~ 1-2 MeV binding. We also 
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emphasize that our calculation neglects second-order double-exchange contri- 
butions, which are expected to be attractive and nearly cutoff-independent. 
Furthermore, the effective mass approximation needs to be improved to treat 
the full momentum dependence of the Hartree-Fock self-energy. 



Our results establish that the 3N force drives saturation for low-momentum 
interactions. However, this does not imply that the 3N contributions are un- 
naturally large. To provide further insight, we list in Table 2 the expectation 
values for the kinetic and potential energies in symmetric nuclear matter. The 
expectation values are obtained by the Feynman-Hellman method, where one 
multiplies the operator of interest in the Hamiltonian by a parameter g, so 
that H9 = (H-0) + gO, and computes (O) = (d(H<>) /dg)\ g=1 . The cutoff de- 
pendence of the expectation value for the kinetic energy (T) reflects the fact 
that correlations in the wave functions are not observable and are reduced 
as the cutoff is lowered. Note that (T) is significantly larger for large-cutoff 
interactions, see e.g., [8]. 



We also show in Fig. 7 the ratio of (Mjn) / '(Mow k) for the triton, alpha particle, 
and nuclear matter at various densities. According to chiral power-counting 
estimates, a natural 3N contribution should scale as (Msn) ~ (Q/A) 3 (Mjn)) 
where Q is a typical momentum scale of the low-energy system. Our results 
for the expectation values show that the 3N contributions are not unnatural; 
for example, the ratio (V 3 n) / (V\ ow k) (A/Zcf) 3 is 0.31, 0.51, 0.53, and 0.35 for 
the second-order calculation at k F = 1.2 fin" 1 and A = 1.6 fm -1 , 1.9 fin -1 , 
2.1 fin -1 , and 2.3 fm -1 respectively. Nevertheless, the role of 3N forces in driv- 
ing saturation means that we will ultimately want to identify a power counting 
for finite density that includes both NN and 3N Hartree-Fock contributions 
at leading order. It is also interesting to note that the deviation from the em- 
pirical saturation point is compatible with a 3N or 4N force contribution of 
order (Q/A) 4 . 



Finally, as was previously found for the binding energies of A = 3,4 nu- 
clei [21], we find that the 3N force may also be perturbative in nuclear mat- 
ter calculations for cutoffs A < 2fm~ 1 . For example, the expectation value of 
(Mowfc) for k-p = 1.2 fm" 1 in the absence of 3N forces are (Mow*;) = —37.5 MeV, 
-37.9 MeV, -38.2 MeV, and -38.5 MeV for A = 1.6 fm" 1 , 1.9 fm" 1 , 2.1 fm" 1 , 
and 2.3 fm -1 respectively. Comparing these with the expectation values in Ta- 
ble 2, we find that 3N effects in the wave functions are small. In contrast, 
the expectation value of the Argonne v±s potential changes by 1.6 MeV after 
resummations when the Urbana IX 3N interaction is included [8]. 
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4 Conclusions and Summary 



In this paper, the nonperturbative nature of inter-nucleon interactions is ex- 
plored by varying the momentum cutoff of a two-nucleon potential. Conven- 
tional force models, which have large cutoffs, are nonperturbative because of 
strong short-range repulsion, the iterated tensor interaction, and the presence 
of bound or nearly-bound states. But for low-momentum interactions with 
cutoffs around 2fm" 1 , the softened potential combined with Pauli blocking 
leads to corrections in nuclear matter in the particle-particle channel that are 
well converged at second order in the potential, suggesting that perturbation 
theory (as in a loop expansion) can be used in place of Brueckner resumma- 
tions. A Weinberg eigenvalue analysis provides quantitative backing to these 
observations. 

We present the first calculations of nuclear matter using the low-momentum 
two-nucleon Vi ow fc with a corresponding 3N force from chiral EFT, which are 
fit in free space and to binding energies of three-nucleon systems with no ad- 
justments when applied at finite density. This combination exhibits nuclear 
binding in the Hartree-Fock approximation and becomes significantly less cut- 
off dependent with the inclusion of the dominant second-order contributions. 
The role of the 3N force is essential to obtain saturation, but the contribu- 
tion to the total potential energy is still compatible with EFT power-counting 
estimates. 

At lower cutoffs around 2fm _1 , the iterated tensor interaction in the two- 
body sector does not play a major role in nuclear saturation, in contrast to 
the conventional wisdom. We emphasize, however, that the relative impor- 
tance of contributions to observables from the tensor force or from three-body 
forces are scale or resolution dependent. Renormalizing the potential to a dif- 
ferent cutoff redistributes these contributions. Our point is that potentials 
with low-momentum cutoffs may be superior for practical nuclear many-body 
calculations. 

We further emphasize that the use of V\ owk is not a replacement for EFT 
but can be used to improve EFT interactions for many-body applications. 
The strategy is to start with a consistently truncated chiral EFT (including 
many-nucleon interactions), using the highest feasible cutoff A to minimize the 
truncation error in matching. Then we evolve the EFT to lower cutoffs (lower 
resolution) using the RG. This running does not, in the present case, involve 
large separations of scale that would justify a truncated derivative expansion, 
so we keep induced operators to all orders to prevent the truncation error from 
growing. 

One could be misled into thinking a higher cutoff implies that the physics is 
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more valid. As stressed by Lepage, one should match a low-energy theory to 
data or an underlying theory with a cutoff close to the scale at which unknown 
physics starts to be resolved [9] . There is no advantage to putting it higher, as 
it can be counterproductive because counterterms (which translates into the 
potential here) must cancel loop contributions from intermediate states that 
are incorrectly represented. 6 But this does not mean that EFT calculations 
cut off at the breakdown scale are optimal for applications to nuclei. Instead, 
one should choose a resolution appropriate to the problem by matching the 
chiral EFT at the breakdown scale (A 3-4 far 1 ) and then evolve downward. 
We find that the decrease in resolution to A ~ 2 fa -1 is particularly important 
for reducing the tensor force. 

The application of chiral EFT in perturbation theory by Lutz et al. [37] and 
Kaiser et al. [38] may be justified in part by our results, although further 
detailed comparisons are needed. We expect, however, that the convergence 
will be improved for smaller cutoffs (A = 3.28 fm _1 in [38]). In contrast to our 
results, which are predictions based on potential fits to scattering observables 
with relatively small residual cutoff dependence of observables, Kaiser et al. 
fine tune the cutoff to simulate the effects of omitted contact terms. Contact 
interactions are included in V\ ovl k and their contributions to nuclear matter 
are important. In addition, 27r-exchange 3N contributions, which are essential 
for saturation in the present approach, are not included in Ref. [38]. It is 
conceivable that a substantial part of our 3N force contributions are captured 
by the explicit Delta isobars included in Fritsch et al. [39], since the low-energy 
constants c 3 and c 4 are to a large extent saturated by the Delta isobar [40]. 
Therefore, a detailed comparision with this work will be very interesting for 
the future. 

The present calculations can be improved in several respects, which is reflected 
in residual cutoff dependence of the energy per particle in nuclear matter. 
First, the current calculations are not complete at second order in many-body 
perturbation theory. We expect that the missing second-order contributions 
are, in fact, more important than omitted higher-order contributions. The 
double-exchange c-terms are expected to give attractive and nearly cutoff- 
independent contributions. Thus, a full second-order calculation is a high pri- 
ority. Other improvements would be a better treatment of self-consistency at 
second order. In parallel, we will use the cutoff dependence of Vi ow fc as a tool 
to develop a power counting at finite density that specifies the appropriate 
truncations and justifies our conclusions. 

With the present three-body force, a better determination of the q coefficients 



6 In other words, the model-dependent cores in the conventional NN potentials 
mandate model-dependent parts in the 3N force to counteract the effects of the 
former. 
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is needed, since C3 and C4 contribute at the MeV level. Of course, it will be 
important to eventually evolve the many-body forces to lower cutoffs at the 
same time as the two-nucleon force is renormalized, as this will generate better 
estimates of higher-order three-body and four-body contributions. In addition, 
the impact of using cutoffs for A < 2/cp needs further investigation. 

Finally, there are particle-hole contributions to consider. These contribute to 
the energy at third order. While we have shown that nuclear matter is pertur- 
bative in the particle-particle channel and that 3N forces naturally provides 
a saturation mechanism, there is still the possibility that a nonperturbative 
summation in the particle-hole channels is required [41]. The lore with con- 
ventional potentials is that these correlations are important for the nuclear 
response but not for the energy of bulk nuclear matter. If this also holds for 
low-momentum potentials, then we can conclude that nuclear matter under 
ordinary conditions can be treated perturbatively in the inter-nucleon forces. 

For applications to astrophysics, however, we may seek higher densities (e.g., 
for applications to the neutron star equation of state). The increase in three- 
body contributions with density above saturation implies that one should only 
do controlled extrapolations to higher density. Controlled extrapolations are 
also needed to reliably describe neutron- or proton-rich nuclei towards the drip 
lines. The use of a variable cutoff and natural estimates of omitted contribu- 
tions is the ideal tool for such investigations. 
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A Pionic three-nucleon force contributions without regulator 

In this Appendix, we provide analytical expressions for the In- and 27r-exchange 
3N force contributions at the Hartree-Fock level, when the regulator /r is ne- 
glected. Analogous 27r-exchange contributions with explicit Delta isobars have 
been evaluated in [39]. This provides an estimate for these 3N contributions 
and a check for the numerical integration with regulator. 
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For the l7r-exchange 3N interaction (Vd), we obtain from Eq. (23) with / R = 1 



Ep ] _ g A c D 3 
A -#^ 288(270' 



32^-72 kiml + Uklm^ 



F 7T ~ ^ "T "Sr 



+ 96 k p m T arctan 



f^)-(36A; F mt + 3m«)l g(l + 



Akf 
ml 



and for the 27r-exchange 3N interaction (V c ), Eq. (26) leads to 



g\ml Cl 9 



A klft 2(2tt)< 



dmJm£ / xdxG s (x,u) 



+ 



+ 



128 



J x 2 dx G 2 v (x , u) 



2 k% 



J xdxG s (x,u) - kp ™ n d m7r I ml J xdxG s (x,u) 



(A.l) 



+ 



(A.2) 



Here u = kp/m n and the auxiliary functions G s (x, u), G v (x,u) and G t (x,u) 
are defined as 



G s (x, u) 



Aux 



(2u 2 — 3) + 4rr(arctan(w + x) + arctan(w — 
1 + [u + xf 



+ {x I -u 1 -l)\og- 

1 + (u — xy 

G v (x,u) = - 2 ( (« 4 - 2u\x 2 - 1) + (* 2 + I) 2 ) log ] + i ( U + % 

x z \ x J i + \ u — xy 

-Aux {x 2 + u 2 + 1)^ 

G,(x, M ) = ^(8 M 2 + 3x 2 )-f (l + « 2 ) 2 
6 2x 
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+ - I ir 2 x 4 + (1 - 3w 2 )(l + w 2 - x 2 ) log 



8 V x 2 



1 + (« - x) 2 

(A.3) 
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